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Structure of monolayers built during adsorption process is strongly re- 
lated to the properties of adsorbed particles. The most important factor 
is their shape. For example, adsorption of elongated molecules on pat- 
terned surfaces may produce certain orientational order inside a covering 
layer. This study, however, focuses on random adsorption of dimers on 
flat, homogeneous surfaces. It has been observed that despite the lack of 
global orientational ordering, adsorbed dimers may form local, orientation- 
ally ordered structures [Tjrj]. Our investigations focus on the dependence 
between domain size distribution and environmental parameters such as 
ionic strength, which affects the range of electrostatic interaction between 
molecules. 
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1. Introduction 

Irreversible adsorption of particles at solid and liquid interfaces is of 
major significance for many fields such as medicine and material sciences as 
well as pharmaceutical and cosmetic industries e.g. 013]. For example ad- 
sorption of some proteins plays crucial role in blood coagulation, fouling of 
contact lenses, and plaque formation. The process is also used in membrane 
filtration units. Additionally, controlled adsorption is fundamental for ef- 
ficient chromatographic separation and purification, gel electrophoresis as 
well as biosensing and immunological assay performance. 
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Adsorption layers formed on flat, homogeneous surface typically exhibit 
random structure without any signs of order. However, a number of latest 
studies on adsorption focuses on deposition on specially prepared, patterned 
surfaces [5-9]. It is due to the possibility of gaining better control over ad- 
sorption process and the properties of resulting mono- and multilayers. In 
such experiments, adsorbed particles take on order of the base. On the other 
hand, anisotropic particles can form orientationally ordered layers even on 
flat, homogeneous surfaces. Ordering can originate from non-uniform align- 
ment of particles in a bulk phase [10] or appear spontaneously in a form of 
ordered domains. The second is the subject of investigation in our study. 
As far as we know, such ordering hasnt been investigated before. Although 
local ordering has been reported in earlier numerical studies on dimer ad- 
sorption [U|2]) not even basic properties of orientationally ordered domains 
have been further investigated. Here, we focus on domain sizes and their 
spatial distribution for different range of repulsive interaction between ad- 
sorbed particles. 



Adsorbed particle is modeled by two identical touching spheres of radius 
ro (see FigfT]). This model was chosen because it is the simlest numerical 



Fig. 1. Model of dimer particle used in Random Sequential Adsorption simulations. 

model of anisotropic particle. Moreover, it has been carefully studied and 
basic properties of the adsorbed layers are well known [HE]. It can be also 
used as a coarse-grained model for some biomolecules [11]. Electrostatic 
interaction between two dimers originates from electrostatic potentials gen- 
erated by individual spheres. In typical substrate this interaction is ap- 
proximated by an exponentially decaying Yukawa potential known also as 
screened electrostatic potential [T2l[13] : 



2. Model 





for L e > 
for L e = 



(1) 



where a and q are particle diameter and charge, respectively; and r is dis- 
tance between individual spheres from different particles. Qq is a constant 
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parameter, eo is a vacuum permittivity and e r characterizes dielectric prop- 
erties of a solvent. Range of the interaction is controlled by parameter 
L e , known also as the Debye screening length. It is a key parameter for 
electrostatic interaction in electrolytes and can be calculated by |14U15j : 



where ks is the Boltzmann constant, T is temperature, e is elementary 
charge, and I is ionic strength of electrolyte solution. 

Adsorption layers were generated using Random Sequential Adsorption 
algorithm introduced by Feder [18]. It is based on repeating the following 
steps: 

I a new virtual dimer is created on a collector. Its position and orientation 
are selected randomly according to uniform probability distribution; 

II. a the virtual particle undergoes overlapping test with its nearest neigh- 



II. b the total electrostatic potential U between the virtual molecule and 
previously adsorbed dimers is calculated using Eq.([T]); 

III if there is no overlap the virtual particle is added to the existing layer 
with a probability P(U) = exp(C//fcsT) (T61TT7]. Otherwise the virtual 
dimer is removed and abandoned. 

The whole procedure is repeated for a specified number of times, which 
guarantees close approach to the jamming limit. The simulation of adsorp- 
tion process was performed for a square collector of 200a-side size and was 
stopped after 6.37T0 8 iterations. For each set of parameters, 20 to 100 inde- 
pendent simulations were performed. Parameters of electrostatic potential 
(PP) were chosen based on [2] to describe typical experimental conditions of 
water solutions. Therefore, the thermal energy ksT at a room temperature 
acts as energy unit and e r = 78. The ro = 9.29 nm provides length scale 
typical to mid sized bio-molecules, and Qq = 6.78 ksTa/e 2 . 



Particle is treated as a domain member when it is close to other domain 
components and its orientation does not differ much from the mean domain 
orientation. However, specific threshold values of distance and angle may 
have significant impact on domain size and shape. Moreover, mean distance 
between particles in a layer depends on electrostatic repulsion range L e . It 




(2) 



bors; 



2.1. Definition of domain 
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means that the same distance threshold value for all L e will prevent domains 
to be classified as such when interaction range is large. On the other hand, 
such "soft" interaction affects the effective anisotropy of a particle which 
changes the tendency for domain formation. Figj2] shows ratio of particles 
classified as domain members in respect to specific values of distance and 
angle difference used in classification process. 



Fig. 2. Ratio of particles in domains as a function of distance and angle difference 
used in cluster classification procedure. L e = 1.0. 

For the purpose of the present study, the concept of domain is defined 
according to Refs. [19i[20j. Particle is domain member when it fulfills both 
of the following conditions: 

I its angle does not differ from mean domain orientation more than 15°; 

II its monomer-to-monomer distance to the closest domain component is 

not larger than 3. 0<r, where a is characteristic size of an interacting 
dimer. Distance is measured between the centres of monomers. 

Parameter a depends on L e . For hard core interaction (L e = 0), a = r$. 
With the growth of L e parameter a increases, to compensate lower coverage 
ratio 9 of a layer: 



For convenience, a values used in this study is collected in Tab J2.1i Using the 
above definition of a cluster, typically, 25% to 45% of dimers are classified 
as domain members. 
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[ro] 


° N 


0.0 


0.547 


1.000 


0.1 


0.510 


1.036 


0.2 


0.485 


1.062 


0.5 


0.435 


1.121 


1.0 


0.372 


1.213 


2.0 


0.277 


1.405 


5.0 


0.146 


1.936 



Table 1. Parameter a ^ for different electrostatic interaction range L e . Values of 
maximal random coverage ratios were taken from [5] . 



3. Results and discussion 

3.1. Hard core interaction 

Adsorption of dimers not interacting electrostatically (L e = 0) provides 
reference level for further investigations. The first important observation is 
that appearance of orientationally ordered dimers clusters is typical for ad- 
sorbed layers generated by RSA algorithm (see Figj3]). The distribution of 




Fig. 3. Fragment of layer build of adsorbed dimers (left) and the same fragment 
but only clusters including at least 5 particles are shown (right). Cluster colour 
describes its mean orientation. L e = 0. 

domain sizes is presented in FigHJ Its exponential decay suggests that there 
is constant probability of classifying particle as a member of orientationally 
ordered cluster. The probability read from plot is equal to 0.532 ± 0.005. 
This value is much higher than 0.167, which is the probability of random 
selection of particle orientation within [—15°, 15°] interval. On the other 
hand, a single cluster borders on more than one dimer, which increases the 
probability of finding nearby parallely aligned particle. Therefore, the above 
results do not answer the question whether domain formation is a merely 
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Fig. 4. Histogram of domain sizes in hard core (L e = 0) interacting dimer layers. 
Dots are measured data whereas solid line is an exponential fit: p(n) = 0.379 ■ 
exp(— 0.631 • n). Insets shows fitted exponentials to domain size distribution for 
partially filled layers. 

statistical effect or rather intrinsic property of elongated particle adsorp- 
tion process. To check it, we analysed cluster statistic before the maximal 
random coverage has been reached. For result comparison purposes, do- 
main classification criteria had to be changed due to lower particle density 
according to ([3]). For example, a = v2 for half of particles inside a layer. 
Results presented in FigJJ] inset show that denser layers are characterized 
by larger probability of parallel alignment. For example, after adsorption of 
only 10% of the particles, the probability of classifying particle as a domain 
component is 0.317 ± 0.006, which is significantly smaller than the value 
for maximally covered layers. Therefore, it can be concluded that Random 
Sequential Adsorption of anisotropic particles leads to formation of orien- 
tationally ordered clusters even in environment with no specific direction. 

3.2. Screened electrostatic interaction 

Anisotropic molecules characterised by soft interactions also form ori- 
entationally ordered clusters. Moreover, the domain size distribution here 
is exponential, too. The probability of attaching particle to cluster de- 
creases with interaction range; however, it has incomprehensible high value 
for L e = 5.0, which requires further study (see Figj5]). Besides the one re- 
sult, obtained exponents are similar to the ones for hard-core interacting 
dimers at earlier stage of adsorption (a > 1.0). It suggests that the main 
factor affecting formation of orientationally ordered clusters is the density 
of particles inside layers. The result is not as entirely expected, as soft in- 
teraction generated by screened electrostatic potential effectively decreases 
anisotropy of dimers, as reported in [2J. On the other hand, soft inter ac- 
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Fig. 5. Fitted exponentials to domain size distribution for electrostatically repulsing 
dimcrs. Insets shows the same exponentials dependence on parameter a for both: 
hard core and electrostatically interacting dimers. 

tion increases distance, at which shape anisotropy can be noticed by other 
particles. As shown, both effects compensate each other, at least for L e < 1. 

4. Conclusions 

Dimer adsorption on homogeneous flat surfaces was investigated numeri- 
cally using Random Sequential Adsorption (RSA) algorithm. Formed layers 
contain orientationally ordered clusters of particles. Their size distribution 
is exponential. The exponent depends mainly on particle density inside a 
layer but not on the character of interaction between particles, at least when 
interaction range is small. 
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